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Abstract 

Small deviations from purely bosonic behavior of trapped atomic Bose-Einstein condensates are 
investigated with the help of the quon algebra, which interpolates between bosonic and fermionic 
statistics. A previously developed formalism is employed to obtain a generalized version of the 
Gross-Pitaeviskii equation. Two extreme situations are considered, the collapse of the condensate 
for attractive forces and the depletion of the amount of condensed atoms with repulsive forces. 
Experimental discrepancies observed in the parameters governing the collapse and the depletion of 
the condensates can be accounted for by universal fittings of the deformation parameter for each 
case. 
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I. INTRODUCTION 



In many physical problems one has to deal with a large number of identical particles that 
are not of a fundamental character but which are known to be composed by a bound system 
of several fermions. Examples include systems of identical atoms, molecules or nuclei. If 
the number of "fundamental" fermions contained in the composite particle is odd, it is a 
fermion, otherwise it is a boson. In many situations the internal structure of the composite 
particle can be ignored and the system as a whole treated as a collection of interacting or 
noninteracting point-like particles. This is the case for instance, in the theories of Bose- 
Einstein condensation (BEC) of trapped bosonic gases Another example is provided 
by electron-hole bosonic states in semiconductors, called excitons 0. In such systems, the 
bosons have internal structure and finite size, i.e., they are composite bosonic particles. The 
rational for neglecting the internal structure of the atoms in atomic BEC is that one is dealing 
with a very dilute system in the trap. The low density regime makes it very improbable that 
the internal structures of the atoms overlap in the trap, since the average distance between 
atoms in typical condensates is several times the size of an atom. For excitons, the situation 
is not as favorable as in atomic BEC and effects of the internal structure of the bosons might 
play an important role 0. 

A departure from the purely bosonic behavior of the atoms in a trap might occur in 
situations that the central density of the condensate grows beyond a critical value. This 
happens, for example, when interatomic attractive forces tend to push the atoms to the 
center of the trap and the zero-point kinetic energy is no longer able to stabilize the system. 
The collapse is expected to occur when the number of atoms exceeds a critical value, N^r, 
leading to an interaction energy larger than the kinetic energy. 

The aim of the present paper is to set up a framework to evaluate the departure from 
the purely bosonic behavior of a BEC of composite particles. A complete theory aimed at 
such a task should include all the possible degrees of freedom for the constituent particles, 
which is in general highly prohibitive from the computational point of view. In the present 
paper, we use a phenomenological approach, making use of the concept of quons 0. Quons 
are particles that are neither bosons, nor fermions, and the quon creation and annihilation 
operators obey a particular algebra that interpolates between Fermi e Bose algebras. The 
quon algebra is in fact a deformation of the Fermi and Bose algebras, and is such that when 
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a parameter (g) runs from —1 and +1, it interpolates between the Fermi and Bose algebras. 

Recently, a systematic way to build a many-body quon state has been discussed 0, ^ 
and a general formula for a normalized many-quon symmetric state [^] has been found. The 
developed formalism allows to apply the quon algebra to describe physical systems presenting 
deviations from an idealized situation, in the same fashion of the more usual Quantum 
Algebras An example is [Q the three dimensional quonic harmonic oscillator which 
describes in a effective way anharmonic effects. The quonic deformation is also effective for 
incorporating correlations in many-body systems. Applications involving the antisymmetric 
subspace (fermion-like particles) can be also considered and are now under way. 

The use of the symmetric subspace, and q close enough to +1, allows to describe in a 
very natural way the departure from a purely bosonic behavior of systems of composite 



bosons. In order to apply such ideas to the BEC in trapped gases, we derive in Section |T| a 
quonic version of the Gross-Pitaeviskii (GP) equation, which we denote qCP. In the limit of 
g = 1, the qGP equation reduces to the usual GP equation, widely used in the hterature 
Initially, in Section we present a brief discussion on the description of a composite boson 
built from two non-identical fermions and provide a connection with the quon algebra. This 
will allow a first estimate of the deformation parameter q in terms of parameters of the 
system [0, 0. In Section |IV| we present numerical results, and compare with observables 
related to trapped bosonic atoms. Conclusions and future perspectives are presented in 
Section |V]. 



II. RELATION BETWEEN COMPOSITE PARTICLES AND QUONS 
A. Composite particles 

Let us consider a composite boson state with quantum number a as a bound state of two 
distinct fermions 

^iio) = $^<fr«i&iio), (1) 

flU 

where ^^'^ is the Fock-space bound-state amplitude, aj^ and are the fermion creation 
operators, and |0) is the no-particle state (vacuum). The quantum number a stands for 
the center of mass momentum, the internal energy, the spin, and other internal degrees of 
freedom of the composite boson. For instance, the composite boson could be the hydrogen 
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atom, where aj^ and 6j, create an electron and a proton respectively. The /i and u stand 
for the space and internal quantum numbers of the constituent fermions. The sum over n 
and u is to be understood as a sum over discrete quantum numbers and an integral over 
continuous variables. 

The fermion creation and annihilation operators satisfy canonical anti-commutation re- 
lations: 

{a^, a^} = {al, 4} = {6^, K} = {6],, bl} = , 

{a^, 4} = 5^,1. {b^, bl} = 5^,^ . (2) 

It is convenient to work with normalized amplitudes such that 

{a\P) = , (3) 

and therefore 

E = ^o.,^ ■ (4) 

flU 

Using the fermion anticommutation relations of Eq. and the Fock-space amplitude 
normalization Eq. (^, one can easily show that the composite boson operators satisfy the 
following commutation relations: 

[A^,A^] = [Al4]=0, (5) 



[Aa, AI] = So^^p - Aaf3 , (6) 

where A^^ is given by 

= Yl f E + E <bib)j . (7) 

One can also easily show the following commutation relations: 

fj,'u' fi'u' 

The composite nature of the bosons is evident from the presence of A^^, which is a sort of 
"deformation" of the canonical boson algebra. The effect of this term becomes unimportant 
in the infinite tight binding limit, i.e. in the limit of point-like bosons. Eq. (§) shows the lack 
of kinematical independence of the microscopic operators and b^, from the macroscopic 
ones A'^s. 
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B. The quon algebra 

The quon algebra is defined by the deformed commutation relation 

A^aI - qA^A^ = 6a,p, (9) 
where q is the deformation parameter of the algebra; A^ annihilates the vacuum 

^a|0)=0. (10) 

The quon algebra interpolates between the Fermi and Bose algebras as the parameter q 
varies from —1 to +1. Polynomials in the creation operators acting on the vacuum form 
a Fock-like space of vectors, i.e., the quonic Fock space. In Ref. it was shown that the 
squared norm of all vectors in the quonic Fock space remains positive definite, provided that 
q ranges from —1 to +1. Note that no commutation relation can be imposed on 
and AajAp. However, as remarked by Greenberg]^, similarly to the case of normal Bose 
commutation relations, no such rule is needed for practical evaluation of expectation values 
of polynomials in A^ and when Eq. (|r^) holds. Such matrix elements can be evaluated 
with the repeated use of Eq. (^ solely; annihilation operators are moved to the right using 
Eq. (P) until they annihilate the vacuum or creation operators are moved to the left using 
the adjoint of Eq. (H) until they annihilate the vacuum. 

Since there is no commutation relations between Aj^,A^^ and Aa.Aj^ all n\ states given by 
the order permutations of the creation operators acting on the vacuum (and carrying different 
quantum numbers) are linearly independent. This is in sharp contrast to the usual bosonic 
or fermionic algebras and important differences between the quonic and usual Fock space 
appears due to this property. As an example, we consider a two quons system occupying 
the arbitrary states 1 and 2 

|i) = 44|o) |2) = 44|o), (11) 

which are not independent states and the overlap matrix is 

(2|1) (2|2)j \qlj 
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One can show that the basis that diagonahzes the overlap matrix is 

\s) = [ (44 + 44) 1 0, ) (13) 

V2v/1 - q 

14 = -j^^={A\Al - 44) I 0), (14) 

where 15) is a symmetric state and \A) is an antisymmetric state. Analogously for a three 
quons system occupying three distinct states 1, 2, 3, the non-orthogonal basis will contain 
the six independent states 4^]^I|0)) "where the indexes i,j,k are determined by the 3! 
permutations of 1, 2, 3. The orthonormal basis may be obtained classifying these states 
according to the irreducible representations of the permutation group of three elements, S^. 

Besides the well known symmetric and antisymmetric states, there are four more exotic 
mixed symmetric states. To shorten the corresponding expressions we adopt the convention 
44^1'!'-') = I'^j^) '^^^ write: states: 

a/1 + 2q^ + 2q + q^ V6 

14 = — ===L_=i= [\tjk) - \jzk) - \tkj) + \jkz) + \ktj) - \kjz)] , (16) 
a/1 + 2q^ — 2q — q-^ vo 

\MS1) = ^ f 4^[Nj'^) - b'^^) + 2K^J) + \jkt) - 2\ktj) - \kjt)l (17) 

a/1 — + g — g'^ V 12 

\MS2) = I 1 \-\i3k) - m + IjH) + \k3i)\ , (18) 

\MSZ) = ] ^ \\i3k) - m - \jkz) + \kjz)] , (19) 

^1 - g"^ - g + g'^ ^ 

\MS4) = ^ =^ [\tjk) + \jtk) - 2\ikj) + \jki) - 2\kij) + \kji)] , (20) 

a/1 — g^ — g + g'^ V 12 

where z, j. A; = 1, 2, 3 . Also, the cases i=j,i = k, j = k and i = j = k are automatically 
included in the expressions of Eqs. (|TB|) to (^) up to a normalization factor which is g- 
independent. Evidently, the above basis states can be built from the well known procedure 
based on the Young tableaux method 0, or through the diagonalization of the overlap 
matrix obtained from all possible order permutations from the state 44^fcl'^^ ' "'■^ fact, 
the g-polynomials which appear in the square roots in Eqs. (15) to (pO]) , correspond to the 
eigenvalues of the overlap matrix and measure the degree of violation of statistics in the 
system. 

Finally, we remark that for the quon algebra one may also define the transition operator 
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Nij, obeying the commutation relations [Q, ^ 

[Nij , aI] = 6,kAl , [Nji , Ak] = -6jkA^ . (21) 

The Nij operator is given by an infinite series expansion in terms of the quonic annihilation 
and creation operators, and its first few terms are given by 

N,^ = A\A, + (1 - Y.^A\A\ - qA\A\){AjA, - qA,A,) + ■■■. (22) 

7 

The number operator Ni is obtained when i = j. 
C. The quon algebra and composite particles 

If one writes q = 1 — x, the deformed commutator, Eq. (P), can be written as 

[A^, A^] = 6af3 - xAiA(3 . (23) 

The similarity of this relation with Eq. is evident. In some sense, the product of fermion 
operators aj^/Ct^ and bjj^bu' weighted by the $'s in Eq. (H) is effectively modelled by the term 
xAlAp in Eq. (|2|). 

We note that the commutation relations of Eq. (|^) between Aa,A/3 and v4j^,yl^ are not 
algebraic relations valid for quons. Since we wish to describe a system of "identical" com- 
posite bosons, we impose that the state vectors of a many-body composite particles system 
are invariant by the permutation of the particle indexes. So we assume that the physical 
subspace which is adequate for the description of a bosonic composite system is composed 
only of totally symmetric states. To build the physical subspace we have to project out from 
the basis of states only totally symmetric states. 

It is possible to show that the most general symmetric state for a system of quons can 
be written as [01: 



Kn,n,..- S) = ^^^^^1^^ (24) 

where is an operator that generates all possible combinations that are symmetric under 
the permutation of any of the creation operators, + rij^ + + ... = N and ||^, 11 



with [A^]! = [A^][A^ — 1]....[2][1] and [0]!=1. Another important result that we are going to 
use next is the following Q]: 

AalriaUfsTi^...; S) = ^J^-^/n^\na - 1, n^n^....; S). (26) 

This last expression allows one to calculate matrix elements between symmetric quonic 
states with any number of quons. Using Eq. (^) one can easily show that, within the physical 
subspace, i.e., the one built only from totally symmetric states, 1^,5"), the commutation 
relations of Eq. are now valid also for the quons 

^1 {AiA\ - ) 1^, ^) = , (27) 

(V',5|(A„A^-A^A„)|^,5) = 0. (28) 
In this way, the analogy between a composite particle and a quon is complete. 



D. Free composite gas 

Now let us consider a system of A^ composite bosons in a large box of volume V at zero 
temperature. If the bosons were ideal point-like particles, the ground state of the system 
would be the one where all bosons condense in the zero momentum state. In the case of 
composite bosons, the closest analog of the ideal gas ground state is 

m = 7=y(^o)''|0) , (29) 

where Aq is the creation operator of a composite boson in its ground state (ground state $) 
and with zero center of mass momentum. Due to the composite nature of the bosons, this 
state incorporates kinematical correlations implied by the Pauli exclusion principle which 
operates on the constituent fermions. Among other effects, the Pauli principle forbids the 
macroscopic occupation of the zero momentum state. The closest analog to the boson 
occupation number in the state of Eq. ( pQ] ) is 

_ {N\AlA,\N) 
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In order to evaluate Eq. ([30|), we consider a spin zero boson and use for the spatial part 
of $ a simple Gaussian form such that the r.m.s radius of the boson is tq. To lowest order 
in the density of the system n = N/V, Nq is given by 

No = N{l- jnri) , (31) 

where 7 is a geometrical factor. The exact value of the geometrical factor 7 depends on the 
form of the internal wave-function of the composite boson. Depending on the exact form 
used for the internal boson wave-function this factor runs from 50 to 250 [g. Note that 
Eq. (0) is valid only for 7nrQ <^ 1, i.e. for extremely low-density systems. For higher-density 
systems, a nonlinear dependence of Nq/N is obtained 0. 



It is apparent from Eq. (31) that in the limit of infinite tight binding, tq 0, one has the 
familiar Bose-Einstein condensation. For finite values of the size of the bosons, the effects of 
the Pauli principle become important and one has a depletion on the amount of condensed 
bosons. Moreover, from Eq. (|3T|) one has that if the size of the bound state is of the order of 
the mean separation of the bosons in the medium, d ~ n~^^^, the depletion is almost total. 
The depletion of the condensation is a direct consequence of the deformation of the boson 
algebra by the term Aajs- 

In the same spirit as in the previous case, we take as the closest analog of the ideal boson 
gas ground state the quon state 

\N) = -^{Alr\0) . (32) 

Note that we are using the same symbols for both the composite and quon algebra annihi- 
lation/creation operators. 

As before, the operator AqAq is the number operator in the zero deformation limit only. 
The effect of the deformation can be evaluated taking the expectation value of AIAq in the 
state |A^), Eq. (|32D . To evaluate the expectation value, we make use of the result given by 
equation Eq. (p5|) 

No = {N\AlAo\N) = [N] = ^-^ . (33) 
Using the above result one obtains to lowest order in x, 

No = i[l - (1 - x)^] ~ A^ M - ^A^x j . (34) 
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Comparing this result with the one of Eq. ( pT)) it is clear that the effect of the deformation, 
for sufficiently low densities, is such that 

Nx ~ . (35) 

That is, the effect of the deformation parameter is proportional to the ratio of the volume 
occupied by the bosons to the volume of the system. Note that for densities not extremely 
low, a nonlinear dependence is expected for the deformation parameter, as remarked 
above. 



III. THE GROSS-PITAEVISKII EQUATION FOR COMPOSITE BOSONS 

We now use the quon algebra formalism in the BEC. We consider a system of composite 
bosons interacting in a spherical harmonic oscillator trap. We assume that the effective 
Hamiltonian describing such system is given by 

H = T + V + Vtrap 

= J2WP)AiAp ^ {c^P\VhS)AiAlA,As + J2Wrap\P)N^p, (36) 

where T, Vtrap and V correspond to the kinetic energy, trap harmonic oscillator potential 
and the interaction among the composite bosons respectively. We take as usual 

T(f) = -|iA„ (37) 

and 

V{x,y)=gSix-y). (38) 

To obtain the equation describing the condensate of composite bosons we assume that the 
lowest energy state is given by 



So we impose the variational principle 



:|0). (39) 



{^|J,S\H\6^p,S) = , (40) 
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where the arbitrary variational symmetrized state has the form 

\6tP,S) = SAlAo\^,S) , (/i^O) (41) 
and the symbol S* is a symmetrizer operator. Eq. (^) is equivalent to the expression 

S\ J2 ( {a\T\P)AlAp + i J2{aP\V\^5)AiAlA,As + {a\V''''^\P)N^/j S) = 0(42) 

From Eq. (^) one obtains the matrix elements 

{ij,S\AiAfs\6^,S) = l|lViV5„o5/5M, (43) 
{ij,S\N^fs\6^,S) = ^6^ohf^{i^,S\AiAlA^As\S^,S) 

= ^^^^ViV(A^-l)(5ao5/3oM50 + 5ao5/3oM<5M)- (44) 

Using these in Eq. P3), one obtains 



[N] ,,,, [N][N-1] 

-jj-Tof, + Vom + j;^ K)oom = 0- (45) 



Therefore, the variational principle leads in coordinate space to the equation 



x) = e(f){x). (46) 



This equation may be interpreted as a generalization of the Gross-Pitaeviskii equation to 
quons and, as mentioned previously will be denoted as qGP. 

IV. THE SOLUTION OF THE QUON GROSS-PITAEVISKII EQUATION AND 
APPLICATIONS 

Many authors have recently presented different methods to solve the usual GP equation. 
Here, we follow a variational approach, in which we expand the ground state wave-function 
in a three-dimensional harmonic oscillator basis with angular momentum I = 0. Good 
convergence has been achieved within the method and we were able to reproduce to very 
good accuracy results obtained in the literature with other methods when we go to the limit 
q=l. 

Initially we consider the case of an attractive interaction, represented by a negative value 
of the scattering length . Of particular interest is the study of the well known collapse, 
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which occurs when the number of atoms in the trap reaches a critical number for which the 
system becomes unstable. In that case, the behavior of the GP solution is characterized by 
the dimensionless constant 

k = (47) 

where bt is the trap oscillator length and a the s-wave scattering length. Several authors (see 



for instance Ref. |jT^ and references therein) have calculated the critical value (fcc) for this 
constant using different methods. The value, for an isotropic trap, was always very close to 
kc = 0.575. Recent measurements indicate a value for kc that is between ~ 10% and 
25% smaller than the theoretical predictions. 

Let us now consider now the qGP equation. According to our discussion in Section (y), 
one may interpret the number of condensed atoms as A^o = [A^]. Eq ( ^61) can be rewritten in 
terms of dimensionless variables in the same fashion as done in the normal GP equation [|I| : 

-Ag + 87rfc(g)|0(x)pj (/)(x) = 2l(f){x). (48) 

Here we have used explicitly the oscillator form for the potential and have written the 
interaction constant in terms of the scattering length a as g = Anh'^a/m. The tilde means 
energy in units of the oscillator energy, 'e = e/hwf, and length is expressed in units of the 
modified oscillator constant , x = x/bt{q), with 

Hq) = bt (^-§y\ (49) 



In addition, 



Obviously the qGP equation of Eq. ( ^81) leads to a critical value kc{q) = 0.575. However, 
in the present case the critical constant to be compared with the experimental value is given 
by 

fc, = 0.575 fMV ' . (51) 



For A'o = [A^] = 2000 and using x independent of A^ and as small as x = 3 x 10"'^ (or 
q = 1 — X = 0.9999), one finds that it is possible to fit the experimental result k^. = 0.522 [jl4|. 
This value of x should be compared to the expression coming from our estimate of the 
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deformation parameter in Eq. (^5l). Using for tq the scattering length |a| ~ 15 A, and using 



the hnear dimension of the trap as given by bt ~ 100 A |]I6| one obtains precisely 

X = 27^ ~ 10-^ (52) 

where for the geometrical factor 7 we used the value 7 = 100, as estimated in Ref. 0. 
This result is of the correct order of magnitude to bring the theoretical value for the critical 
constant within the range of experimental values. 

As a second application, we consider the case of a repulsive interaction. In this case 
the number of atoms in the trap can reach large numbers (N up to ~ 10^), which can 
raise the question of validity for the use of the GP equation for such high densities. In a 
recent calculation |]T3|, effects that go beyond the usual GP (or mean-field) solution were 
also considered and shown to lead to a systematic, small increase in the chemical potential 
of the condensate. In our formalism we may obtain the same behavior by preserving the 
original GP dynamics but relaxing the condition that the particles are true bosons. To see 
how this happens , we now rewrite Eq. (|4^) in a slightly different form: 



Hose + ~ ^1''^^'')'' 

where Hose is the usual three-dimensional harmonic oscillator Hamiltonian. Following our 
previous interpretation that the number of condensed atoms is [A^], we conclude that only 
the term proportional to A^/[A^] — 1 differs from the usual GP equation. However, for q 
close enough to 1, this term is small (<^ 1) even for large values of N, and therefore one 
may treat this term as a perturbation to the non-deformed solution. At this point it is 
worthwhile to note that this term will always increase the energy of the condensate by a 
small amount. Of course, the increase will depends on the value of q (or x) used. In order 
to make a phenomenological estimation of the effect of this perturbation we relate again our 
results of section [II D| for the depletion of the condensate to the one presented in Section III 
of Ref. |jl|. There, an estimation for the depletion is given by: 

-^ = ^VaM^), (54) 

where n(0) is the central (r = 0) density of the condensate. Equating our result, Eq. (p4D, 
with the above and taking in to account that a; <^ 1, we obtain: 
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^/aMO). (55) 



Note that a dependence of the deformation parameter on the number of atoms (or on 
the density) appears very naturally here. As we have discussed in Section [11 D| , an A^- 
independent deformation parameter x is physically reasonable only for low density systems. 
In the present situation of a repulsive interaction, the number of atoms in the trap can be 
many orders of magnitude larger than in the attractive case and it is gratifying to see the 
internal consistency of our formalism, in that it is capable of describing both extremes in a 
very natural way. 

In order to perform numerical calculations, one needs n(0), which can be taken from 
the X = solution of the GP equation and the value for A^, which in turns depends on 
X. However, for our purposes here we may evaluate the order of magnitude of x making 
N = N^'m Eq. (pSj). We also use a/ht = 5 x 10~^. This last value might seem quite big when 



compared to the one used for example in Refs. [0] and |14]. However, recent experimental 
developments [0 allow to change the scattering length by orders of magnitude thanks 
to a technique known as Feshbach resonance. The value chosen here is consistent with, 



for example, the experiment reported in Ref. [|T^]. In Table I we present our numerical 
results for the chemical potential e, for different A'o values. Note that, although our initial 



conditions are different from the ones chosen for the calculations in Ref. |T^, the results 
are qualitatively equivalent, in the sense that we get always a bigger value for the chemical 
potential, as compared to the usual GP solution. 



V. CONCLUSIONS 



We have considered in this work the quon algebra to describe in an effective and phe- 
nomenological way the departure from a purely bosonic behavior of a system of composite 
bosons. The formalism was developed previously and relies on the projection of the whole 
quonic space onto the symmetric subspace. The main idea is to preserve the bosonic be- 
havior and leave to the deformation parameter the description of possible deviations. As 
a specific application we have considered the derivation of the Gross-Pitaeviskii equation 
within the quon algebra, which can be done in a straightforward way using our formalism. 
The interpretation of the modified Gross-Pitaeviskii equation obtained is consistent with 
our initial qualitative considerations. 

As for numerical calculations we have in first place solved our deformed GP equation 
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for the collapse in the case of an attractive force in a trapped bosonic gas. Reinterpreting 
the value of the critical constant for the collapse, we obtained that the deformation of the 
bosonic algebra changes the value of this constant, and the change is in the right direction 
to explain experimental results. Of course, the modification depends on the chosen value for 
the deformation parameter g, which is not easy to estimate at the conditions of the experi- 
ment. We have found that a decrease of order 10~^ compared to the bosonic limit q = 1, and 
independent of the number of atoms in the trap, brings the critical value within the recent 
experimental measurements. Secondly, we have calculated the chemical potential for the 
case of a repulsive force between the trapped atoms. This is a more favorable situation to 
test our model, since recently it became experimentally feasible to create condensates with 
a much larger number of atoms as compared to the attractive case. Using a phenomenolog- 
ical estimation for the deformation parameter, for which a density dependence arises very 
naturally, we are able to calculate the correction in the chemical potential. 

Natural extensions of the present work would include the investigation of the effect of 
the g-deformation on other observables of the BE condensates. Also, there are many other 
interesting systems for which the quon algebra could be useful to describe deviations from 
a purely bosonic behavior of many-body systems of composite bosons. One example, which 
seems particularly interesting in the present context, is the case of excitons which are 
boson-like particles formed by a electron-hole pair in a semiconductor. Work in these direc- 
tions is underway. 
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TABLE I: Chemical potential e - in units of huJt - of the condensate for four values of the number 

of atoms Nq, using the usual GP formalism and the qGP equation. The values of x used in each 
case are also shown and were obtained as explained in the text. 

No ecp eqGP X 

1000 7.31 7.47 1.0 x 10"=^ 

10000 17.84 18.32 1.3 x 10"^ 

100000 44.59 46.69 2.0 x 10"^ 

1000000 111.94 131.16 5.1 x 10"^ 
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